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OPTIMAL DISTRIBUTION OF OPPOSITELY CHARGED PHASES: 
PERFECT SCREENING AND OTHER PROPERTIES 


MARCO BONACINI, HANS KNUPFER, AND MATTHIAS ROGER 

Abstract. We study the minimum energy configuration of a uniform distribution of neg¬ 
ative charge subject to Coulomb repulsive self-interaction and attractive interaction with 
a fixed positively charged domain. After having established existence and uniqueness of 
a minimizing configuration, we prove charge neutrality and the complete screening of the 
Coulomb potential exerted by the positive charge, and we discuss the regularity properties of 
the solution. We also determine, in the variational sense of T-convergence, the limit model 
when the charge density of the negative phase is much higher than the positive one. 


1. Introduction 

In this paper, we investigate ground states of the energy for a system including both 
attractive as well as repulsive Coulomb interactions. The very fundamental nature of such 
nonlocal Coulomb interaction is testified by its ubiquitous presence in nature and by the vast 
number of its occurrences in physical systems. 

We consider the problem of finding the optimal shape taken by a uniform negative distri¬ 
bution of charge interacting with a fixed positively charged region; mathematically, this leads 
to a problem in potential theory (HZIEI]). In our model, minimizing configurations are de¬ 
termined by the interplay between the repulsive self-interaction of the negative phase and the 
attractive interaction between the two oppositely charged regions. We investigate existence of 
global minimizers and their structure, and we obtain charge neutrality and screening as key 
features of our system. In particular, it is noteworthy that the positive phase is completely 
screened by the optimal negative distribution of charge, in the sense that outside the support 
of the two charges the long-range potential exerted by the positive region is canceled by the 
presence of the negative one. 

On the basis of this screening result, we can draw a link to the classical theory of obstacle 
problems PEIEIES]: indeed, the net potential of the optimal configuration can be character¬ 
ized ~ outside the positively charged region and with respect to its own boundary conditions - 
as the solution to an obstacle problem, a fact which in turn entails further regularity properties 
of the minimizer. 

Mathematically, we represent the fixed positively charged domain by a bounded open set 
n"*" C and we are interested in minimizing among configurations 11“ C with finite 
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volume the nonlocal energy 
I{n+,Q-) := 

L L +L L " L L 

Here the first two terms represent the repulsive self-interaction energies of 0"*“ and 17“, re¬ 
spectively, and the third term represents the attractive mutual interaction between fl'*' and 
17“. The present model also arises in the modelling of copolymer-homopolymer blends, see 
the remarks on related models below. 

The first natural question concerns the existence of minimizers for this variational prob¬ 
lem. Since the functional does not include any interfacial penalization, and we just have a 
uniform bound on the charge density, the natural topology for the compactness of minimizing 
sequences is the weak*-topology in While the lower semicontinuity of the functional 

follows from standard arguments in potential theory, a non-trivial issue lies in the fact that 
the limit distribution could include intermediate densities of charge, in the sense that the limit 
function could attain values in the whole interval [0,1]: as a result, the limit configuration 
might not be admissible for our problem. We will however show below that minimizer only 
take values in {0,1}, which allows to bring the negatively charged region and the positively 
charged one as “close” together as possible. 

Our next aim is to identify specific properties of the optimal set. Here we first establish a 
charge neutrality phenomenon: the total negative charge of the optimal configuration equals 
the given positive one, i.e. |f7“| = |f7'’'|. In particular, this shows that configurations with 
nonzero total net charge are unstable in this sense. We also study the case where the total 
negative charge is prescribed: we discuss the minimization of the energy under the additional 
volume constraint |f7“| = A, and analyze the dependence of the solution on the parameter A 
fTheorem 12.2p . proving that a minimizer exists if and only if A < Ihl"*"!. The issue of charge 
neutrality is a central question for systems including interacting positive and negative charges. 
For instance, we refer to the work of Lieb and Simon [26], where charge neutrality is shown 
for minimizers of the Thomas-Fermi energy functional for atomic structures, in the context 
of quantum mechanics. Another related question is whether the maximal negative ionization 
(the number of extra electrons that a neutral atom can bind) remains small: we mention in 
particular the so-called ionization eonjecture, which gives an upper bound on the number of 
electrons that can be bound to an atomic nucleus. For some results in this direction, see e.g. 
|3l[23l[25l|35]. 

A second remarkable property of minimizers is that complete screening is achieved (Theo¬ 
rem [2]3|): the negative charge tends to arrange itself in a layer around the boundary of 17"*“ in 
such a way to cancel the Coulomb potential exerted by the positive charge. Indeed, the net 
potential 

■= / T~r —I dy - / -r-r-—r dy 
Jn+ 4:7 t\x - y\ Jq- ATr\x - y\ 

of the optimal conhguration actually vanishes in the external uncharged region M^\(I7+ U 17“). 
Since it can be proved (under some mild regularity assumptions on the boundary of n"*") that 
(f is strictly positive in the closure of 17"*“, this implies that the positive phase is completely 
surrounded by 17“ and the distance between 17+ and the uncharged space is strictly positive; 
moreover, each connected component of 17“ has to touch the boundary of n"*". Notice that, 
although the minimizer 17“ is in general defined up to a Lebesgue-negligible set, we can always 
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Figure 1 . Sketch of the shape of the minimizer (light-shaded area) corre¬ 
sponding to a given configuration fi"’' (dark-shaded area): the negative charges 
arrange in a layer around the positively charged region. 

select a precise representative, which is in particular an open set, see (12.6p . Such properties 
are established by combining information from the Euler-Lagrange equations with ad hoc 
arguments based on the maximum principle. 

The screening result enables us to draw a connection between the classical theory of the 
obstacle problem and our model, showing that the potential (p oi a minimizing configuration 
is actually a solution of the former. In turn, this allows us to exploit the regularity theory for 
the free boundary of solutions to obstacle-type problems in order to recover further regularity 
properties of the minimizing set (Theorem 12.51) . 

In the last part of the paper, we investigate the regime in which the charge density of 
the negative phase is much higher than the positive one, which is modeled mathematically 
by rescaling the negative charge density by ^ and letting e ^ 0. In this case, we prove T- 
convergence to a limit model where the distribution of the negative charge is described by a 
positive Radon measure (Theorem 16.111 : in turn, we show that the optimal configuration for 
this limit model is attained by a surface distribution of charge on 912"’“ (Proposition 16.2p . 

Related models with Coulomb interaction. Capet and Friesecke investigated in [9] a 
closely related discrete model, where the optimal distribution of N electrons of charge -1 in 
the potential of some fixed positively charged atomic nuclei is determined in the large N limit. 
Under a hard-core radius constraint, which prevents electrons from falling into the nucleus, 
they show via T-convergence that the negative charges tend to uniformly distribute on spheres 
around the atomic nuclei, the number of electrons surrounding each nucleus matching the total 
nuclear charge; in particular, the potential exerted by the nuclei is screened and in the limit 
the monopole moment, the higher multipole moments of each atom, and the interaction energy 
between atoms vanish. Hence our analysis on charge neutrality, screening and on the limit 
surface charge model could also be interpreted as a macroscopic counterpart of the discrete 
analysis developed in [S]. 

Due to the universal nature of Coulomb interaction, we expect that our results could be also 
instrumental in the investigation of more general models, where an interfacial penalization is 
possibly added and the phase D'*' is no longer fixed. Recently, the problem with a single self¬ 
interacting phase of prescribed volume, surrounded by a neutral phase, has been extensively 
studied (see e.g. m El na EH i2oi EZj). 

A minimization problem for two phases with Coulomb interactions and an interfacial energy 
term arises for instance in the modeling of diblock-copolymers. These consist of two subchains 
of different type that repel each other but are chemically bonded, leading to a phase separation 
on a mesoscopic scale. Variational models derived by mean field or density functional theory 
IMl Ea Ea El E] take the form of a nonlocal Cahn-Hilliard energy. A subsequent strong 
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segregation limit results in a nonlocal perimeter problem [33] with a Coulomb-type energy 
contribution. For a mathematical analysis of diblock-copolymer models see for example [U 

132 HQ]. 

In a mixture of diblock-copolymers and homopolymers an additional macroscopic phase 
separation in homopolymers and diblock-copolymers occurs, where now three phases have to 
be distinguished. Choksy and Ren developed a density functional theory a subsequent 
strong segregation reduction leads to an energy of the form 

con^{d{n+UTF)) + cin^{dn+) + C2n^{dn-) + i{n+,n-) ( 1 . 2 ) 

where are constrained to be open sets of finite perimeter, to have disjoint supports 

and equal total volume, and where denotes the Coulomb interaction energy defined 

in dni). In [36] the existence of minimizers has been shown in one space dimension and lower 
and upper bounds on the energy of minimizers have been presented in higher dimensions. 
Furthermore, in m the stability of layered structures has been investigated. 

The model that we analyze in the present paper can be understood as a reduction of ([Ll 
to the case cq = C 2 = 0 and a minimization in Q only, for n"*" given. 

Structure of paper. The paper is organized as follows. The notation and the variational 
setting of the problem are fixed in Section [2l where we also state the main results. A detailed 
discussion of the relaxed model with intermediate densities of charge, instrumental for the 
analysis of the original problem, is performed in Section (3] where the main existence theorems 
are proved. Section |2 contains the proof of the screening property, while in Section [5] the 
relation with the obstacle problem and its consequences are discussed (in particular, we prove 
the regularity of the minimizer). Section [6| is devoted to the analysis of a limit surface-charge 
model. Finally, spherically symmetric configurations are explicitly discussed in the concluding 
Appendix. 

Notation. We denote the ball centered at a point x E with radius p > 0 by Bp{x), 
writing for simplicity Bp for balls centered at the origin. For any measurable set E C we 
denote its Lebesgue measure by |£'| := C^{E). The integral average of an integrable function 
/ over a measurable set E with positive measure is f := f. Sublevel sets of a 

function / are indicated by {/ < a} := {x E : /(x) < a}, and a similar notation is used 
for level sets and superlevel sets. 

2. Setting and main results 

Let n"*" C be a fixed non-empty, bounded and open set. We assume that fl"’' is uniformly, 
positively charged with charge density 1. For any uniformly, negatively charged measurable set 
12“ with finite Lebesgue measure, we consider the corresponding Coulombic energy E{Q~) := 
/(n+, 12“), thus 

=L L +L L 4 ^ " L L ■ 

(2.1) 

Our aim is to find the optimal configuration of the negative charge, under the assumption 
that the two oppositely charged regions do not overlap. We hence consider the minimization 
problem 

min |i?(f2“) ; 12“ C measurable, in"*" n f2“| = 0, |I2“| < oo 


( 2 . 2 ) 
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(notice that we require that has hnite volume in order for the energy (12.ip to be well 
dehned). We also consider the closely related minimization problem where the total negative 
charge is prescribed, which for A > 0 given yields 

min |£'(n“) : Q,~ C measurable, n n~| = 0, |n~| = A | . (2-3) 


The energy (|2.ip can be expressed in different ways. We will usually denote := Xn+ 
and u := Xn-i where Xn± are the characteristic functions of the sets 0^. For given charge 
densities u'^,u, the associated potential ip is dehned as 



u+{y) - u{y) 
47r|x — y\ 


dy. 


(2.4) 


Notice that the potential p solves the elliptic problem 


{ —Ay? = M+ — u, 
limixl^oo W{x)\ = 0 

(see Lemma IQ for the second condition). By classical elliptic regularity, we have p G 
(7i,a(I^3) £qj. gvery a < 1 and p G VFj^’^(M^) for all 1 < p < oo. In addition, p G LP(M^) for 
all p > 3, and Vp G L'?(]R^) for all g > | by [Ml Theorem 4.3, Theorem 10.2]. A standard 
argument, based on integration by parts, shows that the energy of a conhguration can be 
expressed in terms of the associated potential as 


E{n-) = [ |V(pp dx. 


Finally, yet another way to represent the energy is in terms of Sobolev norms: indeed, 

) = ||u^ — "wll^-qRa) ■ 


We now state the main hndings of our analysis. We hrst consider the unconstrained mini¬ 
mization problem (|2.2p . in which the total negative charge is not a priori prescribed, proving 
existence and uniqueness of a minimizing conhguration; interestingly, it turns out that the 
volume of the minimizer matches the volume of the positive charge and the system exhibits 
a charge neutrality phenomenon. 

Theorem 2.1 (The unconstrained problem: Existence and uniqueness). Let 12^ C 6e a 
fixed, non-empty, bounded and open set. Then, the minimum problem admits a unique 

(up to a set of zero Lebesgue measure) solution 12“ C M^. Furthermore, the minimizer satisfies 
the saturation property 

|f2“| = |12+|. 

We also obtain a corresponding result for the case (ESI) when the total negative charge is 
prescribed. 

Theorem 2.2 (The constrained problem: Existence and uniqueness/Nonexistence). Let 
CM.^ be a fixed, non-empty, bounded and open set. Then: 

(i) For every A < |i2^|, there is a unique (up to a set of zero Lebesgue measure) minimizer 
12“ of (ESp. 

(ii) For every A > |I2^|, there is no global minimizer of ESP- 
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We will first prove Theorem 12.21 then Theorem 12. II is an easy consequence of this theorem. 
The main technical difficulty arising when we try to apply the direct method of the Calculus 
of Variations to prove the existence of a minimizer is clearly the following; for a minimizing 
sequence we can pass to a subsequence such that Xq,- u weakly* in L°°(]R^), but 

we can not guarantee that u takes values in {0,1}; in other words, the limit object might 
no longer be a set, with a uniform distribution of charge. This obstacle will be bypassed 
by considering the relaxed problem where we allow for intermediate densities of charge, and 
showing that a minimizer of this auxiliary problem is in fact a minimizer of the original one 
(see Section [3] for the proofs of these results). We remark that a similar strategy was also 
used in [36] for a related one-dimensional model. 

Having established existence of a solution to (12.21) . we now discuss further properties of 
the minimizer. The following theorem, whose proof is given in Section jj] and relies on the 
maximum principle, deals with the screening effect realized by the optimal configuration: the 
associated potential vanishes in the uncharged region. Actually, it turns out that such a prop¬ 
erty (together with the nonnegativity of the potential) uniquely characterizes the minimizer 
(see Remark 14. 2 p . 

Theorem 2.3 (Screening). Assume that O'*" C is a bounded and open set with Lipschitz 
boundary. Let he a representative of the minimizer of (ESD and let (p he the corresponding 
potential, defined in (12.41) . Then p >0 inM.^ and 

(f = 0 almost everywhere in \ (14"'' U Ll~). (2-5) 

After possibly changing Q.~ on a set of Lebesgue measure zero, we have 

Q-= {(^ > 0}\fF. (2.6) 

If satisfies an interior ball condition, then we also have (p > 0 in 14+. 

It is convenient to also introduce a notation for the uncharged region: in view of (12.61) . we 
set 

Qo ■■= m3 \ {(^ > 0} , (2.7) 

which by (12.61) coincides, up to a set of measure zero, with \ (14+ U 14“). Notice that by 
(ESD and (EID we are selecting precise representatives of the sets 14 and 14o, which in general 
are defined up to a set of Lebesgue measure zero, and that with this choice they are open 
sets. 

Based on the screening property and classical maximum principles for subharmonic func¬ 
tions, we establish some further qualitative properties on the shape of the minimizer Ll~. 

Theorem 2.4 (Structure of 14“). Suppose that the assumptions of Theorem \2.3\ hold and let 
14“ be the minimizer of problem (m, given by dM]). Then 14 is open, bounded and the 
following statements hold: 

(i) dist(x,ri^) < 2|14+|^/^ for all X £ Ll~; 

(ii) diaml4“ < (1-|-2\/3) diaml4+; 

(iii) for every connected component V of 14“ we have dV n 914+ / 0; 

(iv) if satisfies an interior ball condition, then dist (rio,914+) > 0. 
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Notice that, as a consequence of Theorem 12.31 and Theorem 12.41 the potential ip of the 
minimizing configuration has compact support. We complete our analysis of the minimum 
problem (|2.2p by discussing some further properties of the minimizer, included the regularity 
of its boundary. This relies heavily on the observation that, as a consequence of Theorem l2.3l 
the potential p associated with a minimizer of (12.2p is in fact a solution to a classical obstacle 
problem. Indeed, as a consequence of the characterization (j2.(i|l . p solves 


I Ay? = xy>o} 

\p>Q . 

It then follows that if solves the obstacle problem 


in 


\r!+, 


( 2 . 8 ) 


mm 


+ 2'ip \ dx : V’ € \ fl"''), V' > 0, if; = p on 


/r3\o^ 
see Proposition 15.11 

The well-established regularity theory for the so-called free boundary of a solution to an 
obstacle problem also yields more information about the regularity of the boundary of flo (for 
a comprehensive account of the available results, see, for instance, the book [32]). 

Theorem 2.5 (Regularity). Under the assumptions of Theorem \2.3l let 12“ he the minimizer 
of problem (1221), let p he the associated potential, and let Hq defined by (12.71) . Then 
p G CjQ’^(M^\n+) and the boundary ofQo has finite Ti.'^-measure locally mM^\f2+. Moreover, 
one has the decomposition dfio = T U S, where V is relatively open in dflo and real analytic, 
while xq & T, if and only if 

idiam ({y? = 0} H Br{xo)) 


lim 

r^0+ 


mine 


= 0 . 


where mindiam(ill) denotes the infimum of the distances between pairs of parallel planes 
enclosing the set E. The Lebesgue density of Qq is 0 at each point ofT,. 

The proof of Theorem 12.51 is given in Section (5] and a more precise characterization of the 
singular points of BQq is given in Proposition 15.21 Notice that the only possible singularities 
allowed in a minimizer are of “cusp-type”, since the set Oq has zero Lebesgue density at such 
points. An example of occurrence of a singular point is presented in Remark lA.31 

In the final section we consider for given u'^ G L^(R^; {0,1}) and for e > 0 the energy 


Eeiu) := 



— u\ 


H-1 


if u G 
else. 


= /k3u< a. 


Here our main result is the Gamma-convergence of Eg to an energy defined on a class of 
positive Radon measures, see Theorem 16.11 Furthermore, in Proposition 16.21 we show that 
minimizers of the limit energy are supported on the boundary and thus describe a surface 
charge distribution. 

Remark 2.6. Although we have restricted our analysis to the physically meaningful case of 
three dimensions with a Newtonian potential, we believe that the methods used in this paper 
can be extended in a straightforward way to obtain the corresponding results in higher space 
dimensions. Indeed, our analysis is based on general tools rather than on the specific three- 
dimensional structure of the problem. Similarly, it should be possible to treat also more 
general Riesz kernels in the energy. 






MARCO BONACINI, HANS KNUPFER, AND MATTHIAS ROGER 


3. Existence and the relaxed problem 

In this section, we give the proofs of Theorems 12.11 and 12.21 In order to overcome the 
difficulties in the proof of the existence of a minimizer pointed out in the discussion above, 
it is convenient to relax the problem by allowing for intermediate densities of charge taking 
values in [0,1], the convex hull of {0,1}. 

In this section, we will always assume that n"*" is an open and bounded set with |n'*'| = m 
for some m > 0. We also recall that := xq,+ is the characteristic function of n"*". We then 
consider, for A > 0, the relaxed minimum problem 


min |f(u) : tt G L^( 


[ 0 , 1 ]), 


u dx = 0, 


u dx < X 


'n+ 


(3.1) 


where 


£{u) := / dx 


and y? is the potential associated to u, defined by 
configurations is given by 


The corresponding class of admissible 


AIa := G [0, Ij) : / u dx = 0, / tt dx < a| . 

1 Jn+ Jr3 ) 

We first note that the potential ip is uniformly bounded and indeed vanishes for |x| 
This is a priori not clear, since u may have unbounded support. 


Lemma 3.1. Assume u G A\. Then the potential ip, defined in (12.4p . satisfies 

for all X G M^, 

, ‘iTV / 

for |x| 


1 /3A\ § / N ^ 1 /3m\ 


dvr/ 

\ip{x)\ 


0 


oo. 


oo. 


(3.2) 

(3.3) 


Proof. For t > 0 let r(t) denote the radius of a ball with volume t, thus -j^r(t)'^ = t. By 


classical rearrangement inequalities [241 Theorem 3.4] we deduce 

u+(y) 


JR3 47r|x - y| 


1 , rimfi 

dy = —— 


4vr|y| 


This shows the upper estimate in ()3.2I) . The lower bound follows similarly. 

Next let e > 0 be given and fix i?;. > 4 such that u < e and fl'*' C . Again by 

rearrangement inequalities we can bound 


f ^(y) j / 

/ -—j-r dy < max max 

dTr\x-y\ X o<w<i_ 

fw<e 

Then for every x with |x| > 2R^ one has 

|n+(y) - u{y)\ 


w{y) 


47r|x — y\ 


dy < 


1 


IB, 


(C 


4vr|?/| 


dy = 


r(c)^ 


l¥^(^)l < [ 


47r|x — y\ 


dy + 


/ 


uiv) 




, m + A r(e)^ 


which shows (13.3p . 


□ 


Existence and uniqueness of a minimizer for the relaxed problem (13.ip follow directly from 
standard arguments. 
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Proposition 3.2 (Minimizer for the relaxed problem). For every A > 0, the relaxed minimum 
problem (IMI) admits a unique solution u\ G .Aa- 


Proof. The existence of a minimizer follows by the Direct Method of the Calculus of Variations 
and standard semicontinuity arguments. Indeed, for a minimizing sequence Un G A\ we 
have that, up to subsequences, Un ^ u weakly* in for some measurable function 

u G which is clearly still an element of the class ^a- To prove semicontinuity, we 

express the total energy as 


£{u) 



u){x){u+ - u){y) 
47r|x — y\ 


dxdy. 


For the self-interaction energy of u we have 

u{x)u{y) 


/ / 

JR3 Jr: 


47r|x — y\ 


dxdy < liminf 


Un{x)Un{y) 
47r|x — y\ 


dxdy 


by classical potential theory (see, for instance, m equation (1.4.5)]). For the mixed term, 
we have 


u+{x)un{y) 


dxdy = / (/?+(y)un(y) dy 
47r|x - y| Jr3 

Jr3 Jr3 , 


dvrjx — y| 


dxdy 


where in passing to the limit we used the fact that the potential (p^ associated to the positive 
phase D'*' is a continuous function vanishing at infinity. This completes the proof of existence. 


Uniqueness of the minimizer follows by convexity of the problem: let ui,U2 G .4a be 
two solutions to the the relaxed minimum problem (|3.ip . and let pi,(p 2 be the associated 
potentials. Setting Ua ■= olu\ + (1 — a)u 2 for a G (0,1), we have Ua G .4a and the associated 
potential is given by ipa = cuLpx + (1 — ol)'P 2 - Hence 

£{pa) = / dx < a / Idx-|-(1 — a) / |V(^ 2 p dx = minT(u) 

Jr3 Jr3 Jr3 Ax 

unless = Vy52. Hence ui = U 2 almost everywhere, and the minimizer is unique. □ 


We now turn our attention to some useful properties of a minimizer of the relaxed problem 
mi, following from first variation arguments. 

Lemma 3.3 (First variation of the relaxed problem). Assume that u is the minimizer of the 
relaxed problem dsn and let ip be the associated potential. Let y be any bounded Lebesgue 
integrable function such that jyj = 0. Then the following properties hold: 

(i) If Jjj 3 y = 0 and there exists <5 > 0 such that suppy C {5 < n < 1 — 5}, then 

/ prj dx = 0 . 

Jr3 

(ii) ^//r3^ < 0 and there exists <5 > 0 such that y > 0 on {u < 5} and y < 0 on 
{u > 1 — d}, then 


/ py dx < 0 . 
R3 
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u < X, r] > 0, and if there exists 5 > 0 with suppry C {?x < 1 — 5} then 

Iff] dx <0. 


j 


Proof. We first prove (i) The function := udzei], for e > 0 sufficiently small, is admissible 
in the relaxed problem (13.ip . Let be such that = rj, so that = ip ± e^p satishes 
—A(/?£ = u'^ — Us- Then by minimality of u we have 

[ iVipl"^ dx < [ iVyjgp dx, 

JR3 JR3 

from which, by letting e ^ 0, we immediately deduce 

0 = / Vip ■ Vfj dx = — ipAif dx = — (pi] dx . 


Let now u, rj satisfy the assumptions in (ii) or (hi) Then the function := u + erj, for e > 0 
sufficiently small, is admissible in the relaxed problem (13.ID . and arguing as before we obtain 

0 < / V(/j • V'i/’ dx = — (pAif dx = — ipi] dx , 

JRS Jr3 Jrs 


which completes the proof. 


□ 


As a consequence of the first order conditions proved in previous lemma, it follows that the 
potential associated to a minimizer is everywhere nonnegative. 


Lemma 3.4 (Nonnegativity of (p). Assume that u is the minimizer of the relaxed problem 
dSU and let ip he the associated potential. Then ip > 0 in M^. 


Proof. For 5 > 0, let x E Es := {u > 5} such that Es has positive Lebesgue density at x. By 
an application of Lemma [3^ii) | with rj ;= —XEsnBr{x)} we then get for every r > 0 


L 


EgHBrix) 


p>{y) dy>0. 


Since ip is continuous, it follows that ip{x) > 0 for all x & Es := {u > <5} such that Es has 
positive Lebesgue density at x. By [29l Corollary 2.14], we hence have ip > 0 a.e. in Es. Since 
(5 > 0 is arbitrary, it follows that ip > 0 a.e. in Eq := {u > 0}. By changing u on a set of 
Lebesgue measure zero, we hence may assume that (p > 0 in Eq. 


By the above calculation, the open set U := {ip < 0} is contained in {u < 0}, and hence 
— Aip > 0 in 17. Since ip vanishes at the boundary of U and at infinity, by the minimum 
principle we conclude that ip must be nonnegative in U, which is a contradiction unless 
U = (h. This shows that (/? > 0 in □ 


The following simple lemma is used in the proofs of Proposition 13.61 and Theorem 12.31 
Lemma 3.5. Let w E L^(M^) n L°°(M^) and let ip he the associated potential, that is 


•= / m-1 ^y- 

7r3 dvr X - y\ 


Then 


[ ip dn^ = [ ^( [ w dx) dr 

JdBu Jr ^ \JBr J 


(3.4) 
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In particular, if suppw C Bji for some R > 0 and L 3 w = X, then 


I 

Jobr 


if dW = \R. 


Proof. Since —Aip = w, we have 


Ir 


^ = L. L =- 4 ^ L 


w dx. 


(3.5) 


(3.6) 


Integrating (|3.6p between R and 00 and recalling that liniij^oo JgBn = 0 by Lemma 13.11 


We next use the first variation formulas to show that minimizers of the relaxed problem (|3.ip 
minimize the absolute value of the total net charge within the set of admissible configurations. 

Proposition 3.6 (Saturation of charges). For every A > 0, the solution u\ of the relaxed 
minimum problem dMD with |n'*'| = m satisfies 


ux dx = min{A, m}. 
Furthermore, for all X > m, we have u\ = Um- 


I 


(3.7) 


Proof. Denote by the potential of u'^—u\ as in (12.41) and choose Rq> d such that D'*" C Brq. 
Arguing by contradiction, we first assume 


yt := / Ux dx < min{A, m}. 

4r3 


(3.8) 


Our argument is based on the fact that screening is not possible under the assumption (13.81) . 
Indeed, we will even show 

\Bji n {(^ > 0} n {ux < 1}| = 00 for all R > 0. (3-9) 

We first note that (j3.8l) yields a contradiction if (|3.9p holds. Indeed, by (13.9p we can choose 
6 > 0, R > Rq such that \{Bb^i \ Br) n {(/j > 0} n {ux < 1 — (5}| > 0. Letting rj := 




(BR+i\B[i)C\{ux<l 


_ 5 }, by Lemma I3.1|fiii) we deduce that 


0 > 


/ (pr]dx = / 

JR3 J(l 


^pdx > 0 , 

(.Bi{+i\Si{)n{tiA<l—<5} 

Ux > min{A, m} and proves (|3.7p for X < m, since 


which is impossible. This shows that 
ux G Ax. 

We next give the argument for (13.91) . By (13.41) . we have for all R > Rq 

Jobr Jr ^ JBr ^ JR ^ 

By integrating this identity from R to 00 , we get (p = 00 . Since p is uniformly bounded. 

Ft 

this implies \B'j^ n {(/J > 0}| = 00 . On the other hand, we have |{ua = 1}| < /rs ux = p, 
which yields (|3.9I) . 

It remains to consider the case X > m and to show that the assumption 


P ■= 


j 

Jr: 


ua dx > m 


(3.10) 








12 


MARCO BONACINI, HANS KNUPFER, AND MATTHIAS ROGER 


yields a contradiction. By Lemma [3.41 we have (^ > 0 in On the other hand, by the proof 
of Lemma 13.51 we get 


dR 



(f 





u\) dx = 


47ri?2 


H — m — 


/]R3\b^ 


u\ dx 


Since by (13.10^ the last term is positive for R sufficiently large, it follows that the mean value 
isSfl of on spheres is a strictly increasing function of the radius, for large radii, vanishing 
in the limit as i? ^ oo. This is clearly in contradiction with the fact that 99 > 0. We conclude 
that J u\ = m for X > m and, in turn, ux = Um by uniqueness of the minimizer. □ 


The previous proposition allows us to draw some conclusions on the dependence of the 
minimal energy in the relaxed problem (13.111 on the parameter A. 


Corollary 3.7 (Minimal energy as function of A). For X > 0 let e(A) := £{ux) = min_ 4 j^ £. 
Then e(A) is eontinuous, strictly decreasing for X E [0,m] and constant for X> m. 

Proof. Since AIa C Ax/ for A < A' the minimal energy e(A) is decreasing. The strict mono¬ 
tonicity of e(A), for X < m, follows from the fact that if e(A) = e(A') for some 0 < A < A' < m, 
then the uniqueness of minimizers would imply that ux = ux ', which is not permitted by (j3.7|) . 
The fact that e(A) is constant for X> m follows also from Proposition 13.61 

To prove that e is continuous, we observe that for X' > X we can use the function ■^ux' E Ax 
as a competitor in the relaxed minimum problem (13.ip . which yields e(A) < £{yUxi). Hence, 
by monotonicity 

e(A) < li^nf£’(^^t6A') = linHnf e(A') < e(A), 

which implies continuity from the right. Similarly, by considering X" < X and comparing with 
E Ax" we obtain continuity from the left. Together this shows that e is continuous. □ 

We next address the proofs of Theorem 12.21 and Theorem 12.11 on existence and uniqueness 
for the constrained and unconstrained minimum problems. 


Proof of Theorem \2.A We divide the proof into two steps. 

Step 1: the case X < m. By Proposition 13.21 there exists a unique minimizer u ol £ va. the 
class of densities Ax- By Proposition 13.61 we have u = A. It therefore remains to show 
that the set {0 < u < 1} has zero Lebesgue-measure. Arguing by contradiction, we assume 
that |{0 < u < 1}| >0. Then there exists 5 > 0 such that the set '■= {5 < u < 1 — 5} has 
positive measure. We set r/(x) := {(p{x) — c)xus{x), where 

c := -I ip{x) dx . (3-11) 

M 

The function rj satisfies the assumptions of Lemma I3..'^n)[ and we deduce 

0 = / (pr] dx = f ip{(p — c) dx f {p — c)^ dx , 

Jk.3 Jus Jus 

which implies p = c almost everywhere in 14^. By (12. 4p and standard elliptic theory, we have 
p E for all p < 00 . From Stampacchia’s Lemma |15l Proposition 3.23], one can 

deduce that V</? = 0 almost everywhere in and then that A.p = 0 almost everywhere in 
U 5 . On the other hand IXp = u in Us, which contradicts our assumption. 
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Step 2: the case X > m. Suppose that there exists a minimizer 0 of (ESI), and let 
u := Xn-- By Proposition 13.61 the unique minimizer ux of the corresponding relaxed problem 
Em is given by ux = Um, in particular we have 

£{um) <£{u). 


Moreover, by the previous step, Um is in fact the characteristic function of a set. For R > 0 
let := Um + (1 — Um)XBii\Bj^-, where R = R{R) > 0 is chosen such that = A, which 

is equivalent to the condition 

dvr q ~ q C 

X = m + —{R^ — R^) — / Um dx . 

3 JBnXBj^ 

We deduce first that R{R) —oo as i? —)• oo and then ^{R^ — R^) = X — m + o(l) as i? —)• oo. 
Since 0+ is bounded, and since Um takes values in {0,1} almost everywhere we deduce that for 
R sufficiently large u^ is the characteristic function of an admissible set for the minimizing 
problem ()2.3p . We claim that £{u^) —>• £{um) < £{u). For R large enough, this yields a 
contradiction to the statement that tt is a minimizer of (j2.3p . To prove the convergence of 
£{u^) observe that an explicit calculation for the self-interaction energy of XBrXBj^ yields 




1 

47 r|x — y\ 


dydx 


c{3R^ + 2R^ - 5R^R^) 
l 5 cR^ 6 ^ + ciR^6^ + + 036^ 


(see (IA.2p in the Appendix), with 6 = R — R = + 0{R~^). This implies that the self¬ 

interaction energy of the annulus vanishes as i? —>■ 00 . By a similar asymptotic analysis one 
shows that the interaction energy of the annulus with the charge distributions u ~^, Um also 
tends to zero as ii —)• 00 . □ 


Proof of Theorem \2. 1[ By Theorem l2.2l there exists a unique minimizer of the constrained 
minimum problem (12.3p for X = m. We claim that is also the unique solution of the 
unconstrained problem (12.21) . Indeed, the conclusion follows immediately from the fact that 
Xn- is the unique solution of the relaxed problem (13.ip for A = m, and by Proposition 13.61 
and Corollary 13.71 □ 


4. Proof of the screening property 

We now turn to the proof of Theorem 12.31 related to the screening of the positive charge. 

Proof of Theorem 12.21 The nonnegativity of cp is proved in Lemma 13.41 We introduce the 
closed sets 

A~ = {x G : \Br{x) n £4“I > 0 for all r > 0} , 

Aq = {x G : \Br{x) \ (14"'' U f4~)| > 0 for all r > 0} . 

Notice that A~ and Aq are independent of the precise representative 14“ of the minimizer: 
indeed, they coincide with the closures of the sets of points of positive Lebesgue density of 
14“ and of \ (14+U14“), respectively. Furthermore, |f4“\A“| = | (R^ \ (14+Uf4“)) \ Ao| = 0. 
We have selected the sets A~ and Aq as the largest possible sets for which we can apply the 
comparison argument in Step 1 below. We now divide the proof of the theorem into four 
steps. 
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Step 1. In this step of the proof, we show that 

sup < inf if . 
Ao 


(4.1) 


Indeed, for (14.11) it is sufficient to prove that (/^(xo) < ^{xi) for every pair of points xq G Aq 
and xi G A~. Define a variation field t] G [0,1]) by 

|B,.(a;o)\(0+un-)| A X G Br{xQ) \ {Q+ U Q ), 

r/(x) := < — |_B^(a;Dnn-| if x G Br{xi) n D , 

0 otherwise. 

Then rj = 0, and, by an application of Lemma ISTHI iil . we hence get 

4 v{y) dy < 4 (p{x) dx. 

J Br(xo)\{n+yjn-) J Br{xi)r\Q.- 

Since (p is continuous we deduce by letting r 0 that p{xo) < p{xi), which completes the 
proof of (|4.1D . 


Step 2. We next show that 

inf 9 ? = 0 . (4.2) 

A- 

Since 99 > 0 this implies by (14.11) that 99 = 0 in Aq) which proves ()2.5p . 

To prove (021) we first consider the case in which A is unbounded. Then there is a 
sequence {xn)n in A~ with |x„| ^ 00 . In view of Lemma ITTI this implies 99(x„) —)■ 0 and 
hence (j4.2p . 

It remains to consider the case when A~ is bounded: let Rq > 0 he such that D'*' U A~ C 
Bf{^. Then (| 3 . 5 |) yields 99 = 0 for every R > Rq, and since 99 is nonnegative we obtain 
that 99 = 0 in \ S_Rq. Next we define the “interior” (in a measure-theoretic sense) of Aq, 
that is the open set 

Aq := {x G : \Br{x) (H (D“ U D“'“)| = 0 for some r > 0} C Aq • 


As ip is harmonic in Aq and 99 = 0 in \ Bf {^^, we deduce that 99 vanishes in the closure of 
the connected component D of Aq that contains \ S_Rq. If dD n A~ 7 ^ 0, we immediately 
obtain (j4.2p . 

It therefore remains to consider the case dD n A~ = 0, and we now show that this case in 
fact never happens. Hence we argue by contradiction assuming that dD n A~ = 0. 

We first deduce that dD C In fact, by the assumption dDr\A~ = 0 any x G dD\d^~^ 

has positive distance to the sets D+ and A~ , which shows that x G Aq and implies that x 0 dD 
as Aq is open. Therefore dD \ SD'*' = 0. 

We claim next that there exist a point xq G cfD"’' with p{xq) = 0 and a ball Bji C D"*" with 
dBfi (H cID'*' = {xo} (we can assume without loss of generality that the ball is centered at the 
origin). In fact (note that at this stage we do not assume an inner sphere condition), choose 
any x* G dD. Without loss of generality we can assume that 

D'*" n Br{x^) = {{y, t) G X M : t > '0(y)} (H Br{x^,) 

for some Lipschitz function 1 /) : —)• M. Since x* 0 A~ by the contradiction assumption, 

after possibly decreasing r we obtain that 99 = 0 on graph('0) (H Br{x^). Choose now any 
xi G D"'' n i?j,(x*) with R := dist (xi, graph(?/;)) < dist (xi, c?i?r(x*)). Then there exists 





OPTIMAL DISTRIBUTION OF OPPOSITELY CHARGED PHASES 


15 



Figure 2. The construction of an interior ball touching 50+ at a point xq, 
used in the second step of the proof of Theorem 12.31 

xo G graph(i/;) n dBji{xi) and we deduce that xq and the ball Br^xi) enjoy the desired 
properties (see Figure [2]). 

Since —Acp = 1 in Bji and since (p is not constant in Bji by Stampacchia’s Lemma m 
Proposition 3.23], the minimum principle shows that y? > 0 in B^. Then the Hopf boundary 
point Lemma m Lemma 3.4] further implies that du(p{xo) < 0 for u = Since (p is of 

class we conclude that (^(xq +tu) < 0 for t > 0 sufficiently small, which is a contradiction 
and completes the proof of claim (14.211 and, in turn, of (|2.5p . 

Step 3. We now prove (|2.6I) by showing that 

1(0" A{^ > 0}) \L1+| = 0. (4.3) 

By (12.51) . we have </?(x) = 0 for almost every x 0 n+ U and hence {(/? > 0} C n+ U 
up to a Lebesgue nullset. It remains to show that the set U = {(^ = 0} D 0“ satishes 
\U\ =0. Indeed, recalling that ip G IFjq’^(M^) for all p < oo, using Stampacchia’s Lemma m 
Proposition 3.23] as in the proof of Theorem 12.21 we obtain = 0 almost everywhere in U 
and then that Aip = 0 almost everywhere in U. Since on the other hand Aip = 1 in [/, this 
implies \U\ =0. The above arguments together yield (|4.3p . 

Step 4- We finally show that 

min (/? > 0 . (4.4) 

under the assumption that 12+ satisfies the interior ball condition. Indeed, if (14.41) does not 
hold, then we have min^+ ip = 0. By the minimum principle and since —Aip = 1 in 12+, there 
is xo G 512+ such that </?(xo) = 0. By the interior ball condition there exists Bfi(xi) C 12+ with 
5Bj:j(xi) (n5I2+ = {xo}. But then we can argue as in Step 2 above and obtain a contradiction 
to the fact that ^p >0. This proves (|4.4I) . □ 

Remark 4.1. A-posteriori we can identify the set Aq used in the previous proof with the set I2o 
dehned in (|2.7I) . To show this we fix the representative (|2.6p for 12“. Since 12+ has Lipschitz 
boundary j512+j = 0 holds and hence 12“ U 12+ = {(^ > 0} U 12+ up to a set of measure zero. 
By Stampacchia’s Lemma m Proposition 3.23] we deduce as above that jl2+ D {(^ = 0}] = 0 
and thus 12“ U 12+ = > 0} up to a set of measure zero. This proves that 

Aq = {x G : ji?r(x) (H {(/? > 0}] = 0 for some r > 0} 
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and since {ip > 0} is open 
Aq = [x^ 


: Br{x) n {(/9 > 0} = 0 for some r > 0} = Oq. 


Remark 4.2. The screening property uniquely characterizes the minimizer, in the following 
sense: there exists a unique set 12“ (up to a set of Lebesgue measure zero) such that the 
corresponding potential is nonnegative and vanishes outside n"*" U 12“. Indeed, assume by 
contradiction that there exist two sets 12^, 12^ disjoint from 12'''. We set ui = Xq-j ^2 = Xq- 
and let ipi and ip 2 be the corresponding potentials characterized by 

{-/^ipi = xn+ - Xn-, 

\lim|xKooV^i(a^) = 0. 

Then one has (fi > 0 and cpi = 0 almost everywhere in \ (12+ U 12^), z = 1,2. Hence, 
—A((/9 i — ip 2 ) = —{ui — U 2 ) and testing this equation by (pi — ip 2 gives 

/ \V{ipi - ip2)\^ dx = - j {ui - U2){ipi - ip2) dx < d, 

2R3 J]R3 

where in the last step we have used the screening property. Therefore ipi — ip 2 is constant and 
since both vanish at infinity we deduce that (pi = (p 2 , which implies | 12 j'A 122 | = 0. 

Using the screening property we now further characterize 12“ and in particular show that 
this set is essentially bounded. 


Proof of Theorem \2.4\ The fact that 12 , defined by (|2.6I) . is open follows by continuity of ip. 
We now turn to the proofs of the other statements. 


Proof of (i). We hrst recall that by the positivity of ip and (|3.2p 

32/3 


0 < ip{x) < 


2 

-m3 


for all X G 


2(47r)2/3 

Lemma 1]. Consider xq G {(/? > 0} \ 12+ and observe 
\ 12+. By the screening property, the 


holds. We now adapt the proof of 
that for every r < dist (xo,12+) we have Br{xo) C 

function w{x) := ip{x) — ^\x — X(ff is harmonic in Br{xo) C {(^ > 0}, and the maximum 
principle yields 

max w > w(xo) = ip(xo) > 0. 

9(i?7.(3;o)n{93>0}) 

Since w{x) < 0 on d{ip > 0}, we obtain 


32/3 

0 < ipixn) < max w < ——-^7777 
- - QBr{xo) - 2(47r)2/3 


2 

m3 


r 

~6 


thus 


< 


3 • 32/3 2 32 

—< -m3 
(47r)2/3 2 


(4.5) 


Letting r dist (xq, 12+) we obtain (i) 


Proof of (a). By using m < |7r(^diaml2+)^ in (j4.5p we easily obtain the estimate in (ii) 


Proof of (Hi). Let U be a connected component of 12 , and assume by contradiction that 
dV n (912+ = 0. Then —Aip < 0 in V and (^ = 0 on 5U, which implies by the maximum 
principle that (/? < 0 in U, which is a contradiction. 
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Proof of (iv). The fact that dflo and Sri'*' 


continuity of ip and Theorem 12.3 


have positive distance is a consequence of the 

□ 


5. Formulation as obstacle problem and regularity of minimizers 

In the following proposition we show that, as a consequence of (12.Sp . the potential p 
associated with a minimizer of (j2.2l) can be characterized as the solution of an obstacle 
problem. 

Proposition 5.1 (Formulation as obstacle problem). Let n"*" he as in Theorem and let 
D := Brq \ where Rq is chosen so that 12+ C Brq- Then the potential p associated with 
the minimizer Ll~ of (|2.2I) is the unique solution to the obstaele problem 

min|y + 2i/)^ dx : G H^{D), > 0, — p G Hq{D)^ . (5.1) 

Proof. Existence and uniqueness of a solution of (|5.ip can be easily established by the direct 
method of the Calculus of Variations, and by strict convexity of the functional. Moreover, 
one can show (see, for instance, m Section 1.3.2] for details) that the solution V’ belongs to 
for every 1 < p < oo, and that it solves the Euler-Lagrange equations 

AiA = X{v»o} in D, 

< Ip > 0 in D, (5.2) 

Ip = p on dD. 

The conditions in (15.21) completely characterize the minimizer of (15.ip : indeed, if £ 

H^{D) were two different solutions of (|5.2p . then for every p G Hq{D) we would have 

Vipi • Vp + r?X{v-i>o}) dx = 0 , 7 = 1,2. 

Using r] = ipi — ip 2 as a test function and subtracting the two resulting equations, we would 
get 

0 = “ ^ 2 ) 1 ^ + i'tpi - V’2)(x{v^i>o} - X{v> 2 > 0 })) - V’ 2 )|^ dx, 

which implies ^pl = V' 2 - Hence, since p itself is a solution to (|5.2p by (12.81) (which, in turn, 
follows from the screening property), we conclude that p is also the minimizer of (15.ip . □ 

We next exploit the connection to the obstacle problem to deduce regularity properties of 
the free boundary. 

Proof of Theorem, \2.5\. By Proposition [5T] the potential p is the solution to the obstacle prob¬ 
lem (15.11) and solves equation (]2.8I) . This problem has been widely investigated and we collect 
below the main results available in the literature, for whose proofs we refer the reader to 
the presentation in the book [32] and to the references contained therein (see, in particular, 

laiZlElEEl). _ 

First of all, by [32l Theorem 2.3] one has that p G \ H+). The free boundary 

r((p) := d{p > 0} has locally finite 72^-measure in \ 12+ by |32l Lemma 3.13]. 
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Moreover, by [321 Theorem 3.22, Theorem 3.23 and Definition 3.24] it follows that r((/?) = 
Treg U Tsing, where Tj-eg is a relatively open subset of r((/3) with analytic regularity l [32l 
Theorem 4.20]), while xq E Tging if and only if 


lim — min 

r^0+ r 


diam ({y? = 0} n Br{xo)) = 0 


1 [321 Proposition 7.1]), from which it also follows that the Lebesgue density of {ip = 0} is 0 
at each point of Tsing. 

Now the properties in the statement follow by observing that clDo C r(</5) and r((/j)\clDo C 
Tsing: indeed, the second inclusion is a consequence of the fact that a regular point xq E Treg 
has a neighborhood in which r((/7) is regular, which implies that xq E d^Q. □ 


In the following proposition we show how points in the regular part T or in the singular part 
S of clDo (regular points and singular points, respectively) can be characterized in terms of 
the blow-up of the potential ip at those points [7|, and a structure result for the singular part 
[8] . A different characterization can be provided in terms of the Ou-Weiss energy functional, 
see [an [38]. 


Proposition 5.2 (Characterization of the singular set of clDo). Under the assumptions of 
Theorem A2.5l let clDo = T U S, where T is the regular part of dQo and S is the singular part 
of dido- The sets T and S can be characterized as follows: for xq E dQo, the corresponding 
rescaled potential ipr,xo ’i-s 


ip{xo + rx) -ip{xo) 

Tr,xo\X) .— ^2 


Then, after extraction of a subsequence, we have ipr,xo Txq in C'[q“(]R^) for every a E (0,1). 
The blow-up function ipxQ has two possible behaviors, independent of the choice of subsequence: 


either ip XQ resembles a half space solution, i.e. 

Txq{x) = \[{x- e )+]^ 
for some unit vector e E or 

Txo{x) = ^X • AxqX 


(half-space solution) 

(5.3) 

(polynomial solution) 

(5.4) 


for some symmetric matrix Axq with TA Axq = 1. Then xq E T i/ and only if (15.3]) holds, 
while Xq E S z/ and only if dsai) holds. 

Moreover, setting for d = 0,1, 2 


:= {xo E S : dimker = d} , 


each set is contained in a countable union of d-dimensional C^-manifolds. Finally, = 0. 


Proof. For a proof of the classification of regular and singular points in terms of the blow-up 
of the potential, see [321 Theorem 3.22 and Theorem 3.23], while for the structure of S see 
[32l Theorem 7.9]. 

We have only to show that the set is actually empty. Indeed, for xq E S one has the 
decay estimate 


\ip{x) - \Axq{x - Xo) • (x - Xo)| < cr{\x - Xol) |x - Xop 
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where cr is a suitable modulus of continuity (see [32l Proposition 7.7]). This property clearly 
implies that if xq E S® we have (/? > 0 in Bj.{xq) \ {xq} for r > 0 small enough, which in turn 
yields xq ^ 5170. □ 


6. A SURFACE CHARGE MODEL 


In this section we discuss the asymptotic limit, in the variational sense of P-convergence 
(see mm), of our charge distribution model when the charge density of one phase is much 
higher than the one of the other: this is achieved by rescaling the negative charge density by 
a factor 7 and by letting e go to zero. In the limit model the admissible configurations are 
described by positive Radon measures supported in \ n"*", with the optimal configuration 
realized by a surface distribution of charge concentrated on . We remark that a similar 
limit model, in the particular case where the fixed domain 17^ is the union of a finite number 
of disjoint balls, was analyzed in [9]. 

Given two positive Radon measures /i, u E A7'''(M^), we introduce the energy 


X(/r,u) := [ 

7r: 


1 


dvrjx — y\ 

and we set := We also define the potential 

1 


d^(x)du(y), 


7rs 


dvrjx — y\ 


dfj-iy), 


( 6 . 1 ) 


and we note that 


(Pf,{x) diy{x) = ip^{y) dfi{y). 
7r3 Jr3 


We will denote by := Xq+-^^ th® measure associated with the uniform charge distribution 
in I7''', and by <y9^+ the associated potential. For A > 0, e > 0, we define the sets 

A\ := |/i E A7^(M^) : supp y \ 17"'', J d^ < a | , 

Ax,e := E AIa : y = uC^, tt : ^ | 


and the functionals on A7(]R^) 

My) ■■= 


+ I{y) if /X E Ax^e 
oo otherwise. 


F{y) := 


-21{y+,y) +I{y) 
oo 


if /u E Ax, 
otherwise. 


Theorem 6.1. Assume that I7'’' C is an open, bounded set with Lipschitz boundary. 
The family of functionals T-converge, as e —>■ 0, to the functional T with respect to 

weak*-convergence in Al(M^). 


Proof We prove the two properties of the definition of P-convergence. 

Liminf inequality. Given /Xg ^ weakly* in Al(M^), we have to show that F{y) < 
liminfg^o .7^e(/^£)- We can assume without loss of generality that liminfe^o l^£(lie) < cc, so 




20 


MARCO BONACINI, HANS KNUPFER, AND MATTHIAS ROGER 


that £ •^x,e and Clearly supp^u C \ O"'', and by lower semicontinuity 

^(M^) < liminfg^owhich implies that ^ G A\. We then need to show that 

-2X(/x+,/r) +l{id) < li^inf (^-2X(/x+,/i£) +X(/X£)^ . 

Since the functional I is lower semicontinuous with respect to weak*-convergence of positive 
measures (see [2T1 equation (1.4.4)]), we immediately have 

^ip) < liminfX(/i£). 

£^•0 


Moreover, the convergence ^ /r and sup^ iJ.s{M.^) < oo imply that 


lim 
€—^0 


[ / d/i£ = / / d/i 

4R3 J]r3 


^) : {|( 7 | > e} is compact for every s > 0}. Hence, since 


for every / G C'q(M^) := {<7 G 
iPfj_+ G C'q(M^) we conclude that 

lim I{fi+, He) = lini / ip + dfie = / d/r = X( 7 i+, ^x), (6.2) 

£^0 Jrs 

which completes the proof of the liminf inequality. 

Limsup inequality. Given a measure pL G A4(M^), we need to construct a recovery sequence 
He ^ H such that limsup^^Q ^^(/re) < J^{h)- We can assume without loss of generality that 
F{h) < 00 , so that H S .4a. 

We first show that without loss of generality we can assume that supp h C C \ . 

Indeed, since 512'*' is Lipschitz, we can define for every <5 > 0 a map G C'°°(]R^; M^) such 
that H"'' CC <I> 5 (n+) and ||<h 5 — ^ 0 as (5 —>■ 0 (the map “pushes” the boundary 

of H'*' in the complement of H'*'). We define the push-forward hs of the measure h by setting 
for every continuous function / 


[ f dns ■= [ f o^s dn- 

JR3 Jr3 


It is not hard to see that hs £ .4a, suppHs CC \ 12+, H5 ^ weakly* in A4(M^) and 
A'ins) —t Tin) as d —)• 0. This shows that it is sufficient to provide a recovery sequence in 
the case supp 71 CC \ 12 +. 


We now reduce to the case of a measure absolutely continuous with respect to the Lebesgue 
measure. Indeed, we define for d > 0 the convolution 

HS := PS* H= / Psi- - y) dniy ), 

2r3 

where ps G C^{Bs), ps > 0, PS = 1 is a sequence of mollifiers. Then hs £ .4a for 6 
sufficiently small (since we are assuming supp/r CC \ 12 +), hs is absolutely continuous 
with respect to the Lebesgue measure and hs ^ P (see [29l Theorem 1.26]). We now show 
that we also have B{hs) ^{p)- Indeed, the convergence of I{h~^iPs) to X[h^iP) can be 
proved exactly as in (| 6 . 2 p : moreover 
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(the first equality is due to the translation invariance of the functional X, while the inequality 
is a consequence of Jensen’s inequality and of the convexity of X), which combined with the 
lower semicontinuity of X leads to = X(/i). This yields T{fis) —t X'(;u). 

Hence, to complete the proof it remains just to provide a recovery sequence in the case of 
a measure ^ G A\ absolutely continuous with respect to the Lebesgue measure and such that 
supp/x CC \ n+. This can be done by a simple truncation argument: denoting by u the 
Lebesgue density of /x, we define /x^ := {u A it is then clear that G Ax^s aiid that 

X'(/x) by the Lebesgue Dominated Convergence Theorem. □ 

In the following proposition we discuss the limit problem, showing that the minimizer is 
obtained by a surface distribution of charge on clD'''. 

Proposition 6.2. Let n G he a solution to the minimum problem 

miniX'(/x) :/X G / d/x = ml. (6.3) 

I Jdn+ ) 

Then /x is the unique minimizer of F over Am- 


Proof. We start by observing that the existence of a minimizer is guaranteed by the direct 
method of the Calculus of Variations. Indeed, given a minimizing sequence {pin)n, by the uni¬ 
form bound /x„((9D'*“) = m we can extract a (not relabeled) subsequence weakly*-converging 
to some positive measure /x supported on SD'*', and such that = m. Moreover, by 

semicontinuity of X with respect to weak*-convergence of positive measures, and by the same 
argument as in (|6.2p with replaced by Hn, we easily obtain that ^ is a minimizer of (j6.3p . 

We now claim that the potential generated by the minimizer ^u, according to (16.11) . 
coincides with the potential outside D+. Fix any point xq G 9D+, let p > 0 and denote 
ap := n Bp{xo)). Consider for e > 0 the measure 


Pe := e'H'^ L (9D’'' n Bp{xo)) + 


m 




which is admissible in problem (|6.3|) for e sufficiently small. Then, by minimality of /x 
0 > -2X(/x+, p) -h X(/x) 2X(p+, pe) - X(pe) 


2eQ;„ 


m 
+ 2s 

- 2e 


/ / 

Jdn+ Ja 

LL 
L 


1 


QQ+ 47r\x - y\ 


dp(x)dp(y) - 


m 


/ / 

Jq+ Ja 


1 


an+ 47r|x - y\ 


dx dy{y) 


dn+nBp{xo) 47r|x y\ 


dx dTL^{y) 


aQ+ Jau+nBp(xo) 47r|x y\ 


dTi"^{x)dy{y) + o(e). 


Dividing by e and letting e ^ O"*" we obtain 

— I {PF - V^m) d/x > — / 
m J9n+ ^ ap Jq^ 


{T^l+ -Ptf) dV? 


V JaQ.+r\Bp{xo) 

and since the right-hand side in the previous inequality converges as p —>■ 0 to (p^+ — p^)(xo), 
we obtain that 
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for every xq G 50"*“. We then conclude that there exists a constant a such that ip := 
ip^+ — (ffj, = a fi-a.e. on 50"'“, and ip < a on 

Observe that, if i?o > 0 denotes a radius such that 0+ C -Brq, by (j3.5|) (which still holds 
in the present setting: see, for instance, m Theorem 6.12]) we have 

/ tp = 0 for every R > Rq . (6.4) 

JoBr 

Now, since ip is superharmonic in \ supp p, ip = a on supp/r and ip vanishes at infinity, by 
the minimum principle we have that 

ip > min{0, a} in \ supp p . (6-5) 


Hence condition (|6.4p excludes the case a > 0. On the other hand, if a < 0 then we would 
have that ip is harmonic in \ 14+, ip < a < 0 on (914+ and ip vanishes at infinity, so that 
(/j < 0 in \ 14+, which is again a contradiction with (16.411 . Thus a = 0 and combining (|6.5I1 
with the fact that ip < a on (914+, we conclude that (/? = 0 on (914+. In turn, this implies 
that (/? = 0 in R^ \ 14+ since ip is harmonic in R^ \ 14^ and vanishes at infinity. We have then 
proved that 


L 


gQ+ 47r|x - y\ 


dp{y) = 


'n+ dvrjx - y\ 


dy for every x G R^ \ 14^ . 


( 6 . 6 ) 


We can now complete the proof of the proposition, showing that p is the minimizer of J- 
over Am- Indeed, for every u G Am we have, using dnsD 

Fiiy) = —2 [ [ - —j -r dx du(y) + f f - —j -r du(x)du(y) 

^ ^ Vrs Jn+ dvrjx - y\ U ^3 dn\x - y\ ^ ^ 

= —2 / [ -—;-r dp(x)du(y) + / [ -—;-r di'(x)du(y) 

As Jan+ 47r|x -y\ 47r|x - y\ ^ ^ 

= / / m-i d(Ai-I^)(x)d(^-z^)(y) - / / —- . dp{x)dp{y) 

Jr3 Jm.3 47r|x - y\ dvrjx - y\ 

= I{p-u)+T{p). 

Using the fact that I{p — i^) > 0, with equality if and only ii p = u (see [211 Theorem 1.15]), 
we obtain the conclusion. □ 

Remark 6.3. The proof of the previous proposition shows, in particular, the following inter¬ 
esting fact: if p solves the minimum problem (|6.3p . then 

I -r-n—I —I x g r^ \ 14 +, 

Jdn+ ^ttJx - y\ Jq+ dvrjx - y\ 

that is, the potential ipf^ generated by p coincides outside of 14+ with the potential ip^+ 
generated by the uniform distribution of charge in 14+. In particular, we again find a complete 
screening property: the potential of /i+ — /r vanishes outside of the support of that measure. 
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Appendix A. Spherically symmetric configurations 


An example in which the shape of the minimize! of (12.2p can be explicitly determined is 
when n"*" has spherical symmetry: if is an annulus, the minimize!’ is given by the union 
of two annuli touching n"*" from the interior and from the exterior respectively (Proposi¬ 
tion [AdJ; in the particular case of a ball, the corresponding minimize! is an annulus around 
(Corollary IA.2D . 

In proving this result we will also compute explicitly the general formula for the energy of 
balls and annuli. We remark that similar expressions were computed in [36| for spherically 
symmetric monolayers and bilayers, where the two phases are adjacent and enclose the same 
volume in any dimensions, by deriving the explicit value of the associated potential. 

In the following, given r 2 > ri > 0 we denote by Cri,r 2 ■= \ the open annulus 

enclosed by the radii ri, r 2 . 


Proposition A.l. Let 12+ = for some R 2 > Ri > 0, and let R^ > 1 be the unique 

solution of2{Rl — 1) — (2(i2^ — 1))^/^ = 0. Let Ll~ be the minimizer of (12.211 . Then: 

(i) If ^ < R^ thenll~ = UC'K 2 ,r 2 ; where ri G (0,i?i) andr 2 > R 2 are determined 

by 

irf- rl + 2{Rl - R\) = 0, 

\rl-rl + 2{Rl-Rl)=Q. 

(ii) If ^ > R* then Il~ = Br^ U CR^^r with r = -^2(i?2 “ Bf)- 
Proof We divide the proof into two steps. 

Step 1. We start by computing the self-interaction energy of an annulus Cr^,r 2 ™ spherical 
coordinates: 


(A.l) 


Cr-^,r2 


47r|a: — 


1 rr2 r^TT ^TT / r 

-. dxdy= / dp d9 dy. / 

■ — y\ Jri Jo Jo \JCr^, 


p sm ip 


4:TT\x{p,e,ip) - y\ 


dy 


where x{p,6,ip) denotes the point in Cr^^r 2 whose spherical coordinates are {p,6,ip). Since 
by rotation invariance the inner integral depends in fact only on p, we can compute it for 
X = pes, obtaining 


LI 


1 


r dvrla: — y\ 

'—'7’^ ,7*2 ' ^ ' 

cr2 / i'r2 rSTT 


dxdy = 




dy 


Cr,.ro \P^3-y\ 


dp 


P 


I'Z'K PT 

Jo Jo 


r sm 


= 27r 


'ri JO JO (r^-I- _ 2 rpcos (/?)^/2 
rr2 rr2 

/ / pr (r-|-p — |r — pI) dpdr . 

J r^ J r\ 


dpd^dr jdp 


'ri Jri 

From the explicit computation of the last integral we conclude that 

1 , . dvr 


LI 


r dvrlx — y\ 

V_.'7'-| f'n I ty I 


dxdy = ^ ^3rf -|- 2r| — 5rfr|^ . 


(A.2) 


■'ri,r2 ^ ^ri,r2 

By similar computations we obtain the interaction energy of two disjoint annuli: for 72 > 
Fi > i ?2 > ^1 > 0 we have 
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Figure 3. Example of occurrence of a singularity at the origin in the bound¬ 
ary of a minimizer. 

Step 2. We now turn to the proof of the statement. Notice that by scaling we can assume 
without loss of generality that the inner radius of the annulus is equal to 1, and in particular 
we can consider n"*" = Ci^r for some R > 1. By uniqueness, 0“ is invariant under rotations, 
so that it consists of a union of annuli. By Theorem 12.41 14“ is the union of an annulus, 
touching 14+ from the exterior, and an annulus or a ball touching 14+ from the interior. 

The hrst possibility is that the connected component of 14“ internal to 14+ is not the unit 
ball; the minimizer is then the union of two annuli, 14“ = U Cji^r 2 - The conditions (lA.ll) 
are induced by the volume constraint and by choosing ri > 0, r 2 > 1 such that the energy, 
computed in the previous step, has a local minimum (or, equivalently that ip vanishes on dBr^ 
and dBr^)- Since by (|A.1I) 

2{R^-l) = rl-rl < = (2(i?3_i))§ 

2 

we conclude that this possibility can only occur for 2{R? — 1) — {2{R? — 1 )) ^ < Q, which is 
equivalent to R < R^. 

The second possibility is that the connected component of 14“ internal to 14+ is the unit 
ball, with an external annulus where r is determined by the volume constraint. After 

some calculations we obtain in this case that 

</p(0) = (442 _i)_ 1(2(443 _!))§. 

Since for any minimizer (^ > 0 holds by Lemma 13.41 we obtain that the second possibility can 

2 

only occur if 2{R^ — 1) — (2(44^ — 1)) ® > 0, which is equivalent to 44 > 44* (in case of equality 
we just have ri = 0 in the first case, that is the inner annulus degenerates to the unit-ball). 
Since we always are in one of the two cases the claim follows. □ 

Corollary A.2 (The case of a ball). Let 14+ = Bji for some 44 > 0. Then the minimizer of 
(12.211 is the annulus 14“ = 

Proof. The result follows immediately from Proposition lA.ll bv letting 44i —)• 0. In this case 
the argument is actually simpler and we can give a short independent proof: by Theorem l2.ll 
the minimizer is unique up to a set of vanishing Lebesgue measure. Since 14+ is invariant with 
respect to rotations, it hence follows that also 14“ is rotationally invariant and hence consists 
of an union of annuli. The annulus 14“ = is the only set with this property such 

that Theorem 12.41 holds for every connected component V of 14“ and satisfying the charge 
neutrality condition |14“| = |14+|. □ 
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Remark A.3. We can construct a simple example of occurrence of a singularity in the boundary 
of a minimizer as follows. Let := Bi{\/2ei) U Bi{—\/2ei) be the union of two disjoint 
unit balls centered at points at distance 2^/2. The corresponding minimizing configuration is 
then the union of the two annuli which minimize separately the energy for the two connected 
components of 0+, i.e. Q.~ = Indeed, by linearity the resulting 

potential is nonnegative and vanishes outside n+Ufl”, and by uniqueness of the configuration 
having the screening property (see Remark 14.21) we conclude that kl~ is the minimizer. Notice 
that the two annuli touch at the origin, which is thus a singular point (see Figure [3|). 
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